Optical conductivity of graphene in the presence of random lattice deformations 
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We study the influence of lattice deformations on the optical conductivity of a two-dimensional 
electron gas. Lattice deformations are taken into account by introducing a non-abelian gauge field 
into the Eucledian action of two-dimensional Dirac electrons. This is in analogy to the introduction 
of the gravitation in the four-dimensional quantum field theory. We examine the effect of these 
deformations on the averaged optical conductivity. Within the perturbative theory up to second 
order we show that corrections of the conductivity due to the deformations cancel each other exactly. 
We argue that these corrections vanish to any order in perturbative expansion. 

PACS numbers: 73.22.Pr, 72.80.Vp 



I. INTRODUCTION 

Graphene, a two-dimensional sheet of carbon atoms 
forming a honeycomb lattice, has outstanding electronic 
properties^— . This is due to the fact that there are two 
bands that touch each other at two Dirac nodes. More- 
over, the low-energy quasiparticles of undoped graphene 
experience a linear dispersion around two Dirac nodes. 
Transport properties, characterized by the longitudinal 
conductivity at the Dirac nodes, are quite robust and do 
not vary much from sample to sample. Exactly at the 
Dirac point a minimal conductivity has been observed 
in a number of experiments^ 2 -^. There are two impor- 
tant questions regarding this minimal conductivity: (I) 
is the value of the minimal conductivity "universal" (i.e. 
independent of additional modifications of the graphene 
sheet such as ripples or impurities) and (II) what is its 
actual value in units of e 2 /hi A discrepancy between the 
calculated conductivity of Dirac fermions and the exper- 
imentally observed minimal conductivity of graphene by 
a factor of roughly 1/tt has been the subject of a substan- 
tial number of publications. The central idea is that ei- 
ther disorder— ~— or electron-electron interaction^— may 
affect the value of the minimal conductivity. Moreover, 
the value of minimal conductivity at low temperatures 
depends on the order of varies limits (e.g. frequency 
lu — s- and temperature T — > 0) and is related to the 
scaling property a min {u,T) = a min (oj /T)lk Below we 
will employ the zero-temperature formalism which sug- 
gests T — > and oj 0. This yields for the DC limit of 
the AC conductivity the value tt/2£&^. 

An additional problem in terms of disorder is that it 
is not clear what role is played by different types of dis- 
order. Since disorder, depending on its type, may break 
different internal symmetries of the Dirac Hamiltonian, 
a classification according of the different types is crucial. 
On the other hand, the origin of disorder in graphene 
can be different. Besides impurities inside the graphene 
sheet and in the substrate, the deformation of the lat- 
tice (e.g. ripples) might be the main source of disor- 
der— In general, it is believed that surface corruga- 
tion o 11 ' 17 ' 18 may influence electronic transport properties 



of graphene. It is crucial to notice that lattice defor- 
mations do not break the chiral symmetry at the Dirac 
point, in contrast to potential disorder or a random gap 
caused by a random deposition of hydrogen^. There- 
fore, it is expected that this type of disorder has a rather 
weak effect on transport properties 2 -. This is supported 
by calculations, where the lattice deformations are ap- 
proximated by an uncorrelated random vector potential 
in the Dirac Hamiltonian 21 . This type of disorder has no 
effect on the minimal conductivity-. More recently, how- 
ever, a more general theory of lattice deformations with 
long-range correlations revealed a dramatic increase of 
the minimal conductivity for weak disorder—. In this 
paper we will study a similar model by an alternative 
approach to check whether or not this dramatic increase 
of the minimal conductivity can be reproduced. 

First we consider the deformation of the graphene sheet 
in three dimensions and show that in the continuum limit 
the dynamics of the electrons on the two-dimensional 
surface is defined by the so-called induced Dirac ac- 
tion presented in Ref. [22[ . In our approach the inter- 
nal deformations of the graphene sheet and the defor- 
mations perpendicular to the sheet direction are unified 
into one schema, while in the approach developed in pa- 
per s 1 ^ 21 ' 23 " — there are separate internal 2D gravity and 
additional non-abelian gauge fields. The deformations of 
the sheet in three dimensions by local SO(3) rotations 
of the basic vectors in our approach carry the degrees of 
freedom of the additional gauge field. 

Then we develop a replica-trick based field theory to 
take the random character of surfaces into account and 
to calculate the average optical conductivity by a pertur- 
bative expansion. Our main result is that the random 
lattice deformations do not affect the robust character of 
minimal conductivity, contrary to the result presented in 
Ref. fnl. 



II. THE MODEL 

We depart from a model of hopping fermions on the 
regular 2D honeycomb lattice. Honeycomb lattice has 



FIG. 1: Hopping vectors on regular honeycomb lattice 



natural partition into two triangular sub-lattices and we 
mark electronic fields associated with sites of the sub- 
lattices as (i/jn,ay ^rJ,a)) ot = 1,2. The action of electrons 
hopping on a line with the lattice spacing |e| reads 

t,n 

but when fermions change hopping direction in two di- 
mensional space they fields should also be rotated by a 
corresponding angle (Fig. [TJ . On the honeycomb lattice 
(Fig. [5J we have 

S% *] = t £ (*t,nd t % tH + *t,fl72*U-s!) , (1) 

t,n,i 




i= 1,2,3, 



where 7q , 7j , J = 1,2 are Dirac matrices which are re- 
lated to usual Pauli matrices via 7^ = (T 2 , 7 2 = a 1 , and 
7 = ct 3 and fields 

= = e glS *n, (2a) 

%+s 2 = e^°*n+e 2 = e^T-e^-^fl, (2b) 
*' s+(?3 = e"^T°* ?1+e - 3 =e-^°e^-% s , (2c) 

are rotated by ±47r/3 and translated by 132,3 spinor repre- 
sentations of the rotation group 50(3). In the paper by 
SemenofS^ was shown that the spectrum of low-energy 
excitations of the hopping fermions on honeycomb lattice 
(corresponding to the continuum limit of the model) co- 
incides with the spectrum of Dirac fermions in 3D space. 
Below we will show that the continuum limit of the ac- 
tion of fermions hopping on honeycomb lattice Eq. ([1} is 
defined by the Dirac action in three dimensional coordi- 
nate space. This will allow to construct the continuum 
limit of the generalized hopping model on the randomly 
deformed lattice. 



FIG. 2: (Color online) Regular honeycomb lattice. 

In order to find a continuum limit of the action Eq. (TT]) 

one expands translational operators e £i ' 9 ~ 1 + • d and 
substitute Eqs. (|2a)) for >F's into the action Eq. ([1]). Then 
after some simple algebra one will obtain: 

S[% 9] = iJ2 *t,a U + \l% ■ fd - t]\ (3) 

-> i J d 2 t;dt * U + Iffr ■ fd - t]j 9, (4) 

where we introduced orthonormalized vectors ^ = (e*2 — 
e3)/\/3 and ft 2 — e%. In the line ((U) we have rescaled 
the fields and coordinates as 9 -t 2/3*, £ 3£/2, 
t — > t and fi t — > fj, i . It is clear from Eq. (£[]) that vec- 
tors pi a — /Lt„(^)ei, a — 1,2, with ej, i = 1,2 rep- 
resenting an orthonormal basis in the flat space play 
role of tetrads (vielbein) in a 2D plane with arbitrary 
coordinates £j, i — 1,2. Indeed, consider deforma- 
tion of the honeycomb lattice (cf. Fig. [3]) and attach 
to the sites a new coordinates Then the vectors 
fj^, a — 1,2, i = 1,2 will be connected with the same 
vectors in the old coordinate £j, i — 1,2 via 

/4(l) = /4(f). (5) 

We regard now the vectors /if as vielbeins in a 2D plane 
which obey the orthogonality relation fJtf(J> a ,j — §ij and 
define the metric fif/J," = g%j. After integration by parts 
in Eq. ((3| and using the relation fiiftj = g^ + -^=ejj7 
with [ii — 7 a /i° and g = det[<7,j] one will obtain 

s[% 9] = 1 J d 2 tdt * (a t + 7 v„ [d 3 - \i<Fi\) % (6) 

where Tj — e a6Ma^ j^k,b is a standard spinor connec- 
tion corresponding to the vielbein \i? a and denotes a 
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FIG. 3: (Color online) Random honeycomb lattice. 



covariant derivative. For a scalar function / it reduces to 
a usual partial derivative: V,/ = d i f, while for a vector 
valued function f it is = djj + T^/fc, where T*j 

represent Christoffel symbols. 

Let us now consider deformations of the honeycomb 
lattice in a three dimensional spac e 22 ' 26 . This means that 
two j a , a — 1,2 matrices in a tangent plane become 7 
matrices in a SO(3) rotated plane which is tangent to 
curved surface at the point 

= u&rH.u^). (?) 

As it is shown in Ref. [22| local rotations by U(£) produce 
a 2D surface embedded into 3D Euclidean space if 

U- l dJJ = - A {x a d ll x a + nd^n), (8) 

where x a — ^dfj_x a j a and h — n Q 7 Q , a = 0,1,2 be- 
ing the tangent and normal to the surface 3D vectors at 
the point £ respectively. This will occur since Eq. © 
will fulfill the Gauss-Codazzi equation s 26 ' 27 , which rep- 
resents the necessary conditions for the surface x a to be 
embedded into 3D Euclidean space. 

Then we should rotate also fermionic fields by the same 
matrices ^ —> [/<F, after which the action becomes 

s = ij d\dt w- 1 U + 7 <X[^ - i 7o r,j) UV, (9) 

where T M = iTr(n£/~ d^U) is the spinor connection on 
the surface x a . By use of Eq. (|7|) this expression can be 
simplified essentially (see details in Refs. f22jj26|) acquir- 
ing the form 

S[%*]=i J ^dtffe + ^fi-^jf, (10) 

where 



is the metric on the surface x a (£ ) induced by its embed- 
ding into 3D Euclidean space, g denotes its determinant 
and 

Y = d^x • 7 = d^ l x a la (12) 

represent the induced Dirac matrices^ 2 -. In a flat space, 
i.e. for x(£) = xo + e^ 1 ^ the induced metric reduces to a 
usual diagonal matrix. One can call action Eq. (|T0j) the 
induced Dirac action since the matrices 7 M are induced by 
embedding. The expression in Eq. (|10[) is a generalization 
of 2D action Eq. (gj) to 3D space. 

III. EFFECTIVE ACTION FOR SMALL 
CORRUGATIONS 

Performing integration by parts in Eq. (I10[) we arrive 

at 

5[*,*] = i J d 2 ^*(9 t + V5r[^ + r M ])vI/.(l3) 
Here, the quantity 

r„ = \rv,>% (1 4 ) 

plays the role of an induced spinor connection, where 
denotes the operator of covariant differentiation and is 
defined as 2 ^ 6 . 

v M (--0 = -^(Vs---)- 

Let us derive the asymptotic action for small corruga- 
tions of the graphene sheet. In this case the surface x 
can be asymptotically represented as 

xte 1 ,£ 2 )»x + e%+x / (£ 1 ,&), (15) 

Plugging Eq. (|T5|) into Eq. (|lip we obtain the asymptotics 
of the metric tensor: 

9ufi ~ + e u^ + <>, (16) 

where 

£»n = K ■ 9 M x'. (17) 

Thus the metric tensor is in general neither diagonal nor 
symmetric. Its determinant is found using common rela- 
tions 

S « 1 + 2e u + 2e 22 = 1 + 2e w , (18) 

and correspondingly its square root: 

+ W (19) 

Using Eqs. (fl2|) . (fTSl ) and (fTTjf we arrive at the effective 
action for small fluctuations e: 

S [% e] w ij d 2 £dt * ( 7o 5 t + [1 + ejffd^ * 

+ l -j d 2 Cdt Vyd^*, (20) 
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with induced 7— matrices 7^ = 6^7°. For further pur- 
poses we associate the spatial fluctuations with a bosonic 
field 

<W = A(£t), 

and its gradient with a static vector-disorder like term: 

d ll e vv =d ll \{lt)=B ll (it). 
Hence the action formally becomes 



5 [*,*;A,B] 



e^ctt* (70^+7^)* 



t / d 2 £dt A*7^a„* 



The action derived this way reproduces the ansatz action 
considered in Ref. fl7j . 

Below we consider topologic defects in the flat space. 
Technically that means that we replace zweibeins by 
a unity-matrix. We are ultimately interested in the effect 
of this sort of the disorder on the optical conductivity. In 
order to perform such calculations we have to make some 
suggestions regarding the correlators of the introduced 
quantities. One usually requires the vector disorder fields 
to be gaussian correlated, i.e. 



(22) 



(B^,t)B v (i',t')) = g^jmmz-?), (23) 

which guarantees that the vector associated with the ran- 
dom disorder is static. In Fourier-space these expressions 
read with the short-hand Q = (q Q ,q): 



(BJQ)) = 0. 



(24) 



(BJQ)B V (Q')) = g 2 (2Tr)Xj(q )5(Q + Q'). (25) 



From Eqs. ([23]) and ([25]) we are lead to the correlators of 
the scalars A, since we have an exact relationship 

{B^{Q)B V {Q')) = i\q' u (A(Q)A(Q')), (26) 

which leads to 

(A(Q)A(Q')) = - 1[ ^(2K)H(q )6(Q + Q'), (27) 
q z + [i A 

where we have introduced an infrared cutoff [i of the order 
of the inverse lattice spacing in order to avoid long wave- 
length divergences. Inverse Fourier transform yields for 
the (AA)— correlator 

(A(^)A(eV)) - g^m^iog^d-e)] . (28) 

Furthermore we will always assume 

(A(£t)) = (A(Q)>=0. (29) 



The Fourier transform of Eq. (|20|) expressed in terms of 
the scalar fields A only reads 



A] 



*q(9o7o + 7-9)*q 



A P ^ P+Q T(P + Q,P,Q)^ Q , 



(30) 



Q p 



where 7 • q = 7 M <z M - The two- particle vertex is obtained 
from Eq. (|30[) in limit A — > by performing second or- 
der functional derivative with respect to the Grassmann 
fields: 



(21) where 



T(P,Q) = (2nfS(P-Q)Go 1 (Q), (31) 



Go HQ)= 9o7o+7-9, (32) 



represents the inverse free propagator and correspond- 
ingly 



G (Q) = 



%io + 7 • g 
ql + q 2 



(33) 



the free Dirac propagator. For the three-particle vertex 
function T(K^,K A ,K W ) follows from Eq. {30) 



x ^7-(fe A + 2fc*). 



(34) 



Furthermore we have to augment Eq. (|30p by the in- 
teraction between fermions and the radiation field 



S opt [%*,A] = -J J A p ^ p+q1o * q , (35) 
P Q 

which suggests the presence of an electric field applied 
to the graphene sheet. Interaction Eq. (l3"5"j) gives rise to 
the optical conductivity due to polarization of the charge 
carriers. The corresponding bare vertex is defined as 

T Q (K 9 ,K A ,K 9 ) = (2tt) 3 S(K^ - K A - Jf s ) 7o . (36) 

The full action acquires the form 

5[* > *,A > A] = 5[* J * J A]+5 qpt [* > *,A]. (37) 

The electron-gauge boson interaction renormalizes elec- 
tronic spectrum and therefore should have an effect on 
the response to the radiation field. 



IV. OPTICAL CONDUCTIVITY OF 
GRAPHENE 

We consider first ideal graphene. The corresponding 
Eucledian action is obtained from Eq. (1371) if we assume 
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P+K 




FIG. 4: Bare polarization bubble 

P+K 




FIG. 5: Dressed polarization bubble 



A = 0: 



*q(9o7o+7'9)*q 



J J A{P-Q)^> pl ^ Q , 



(38) 



The optical conductivity of 2D Dirac electron gas can 
be calculated from the electronic polarization^^ 2 - 



7o = 4fc °l i 2 1 o^ nW ' 



(39) 



where Tl(K) denotes the irreducible polarization. Fac- 
tor 4 in front of this expression arises from the taking 
both spin and valley degeneracy into account. To the 
leading order it is given by the diagram shown in Fig. 2) 
Algebraically we have for the polarization bubble: 



n(K) 



Tr{ 7o Go(P)7 G (^ + P)}. 



(40) 



with the bare Dirac propagators G (Q) defined in 
Eq. (|3"3"1) . We give some details of the calculation in the 
Appendix. The irreducible polarization is obtained as 



n(K) 



i 



16 

and the optical conductivity in Si-units is 

D 2 



4 h 



tt e 
2~h 



(41) 



(42) 



In what follows we calculate corrections of the conduc- 
tivity in Eq. (|42p due to lattice deformations described 



in Sec. Mil To the leading order in momenta the renor- 
malized inverse fermionic propagator can be written as 

G-\Q) = GoHQ) - E(Q) « Zo\l + ^7 ■ 9, (43) 
with renormalization factors 

d 



J o.i 



1 



7o, M 
d 



-S(Q) 



Q=0 



(44) 



Q=0 



where £(Q) denotes the fermionic self-energy. On an- 
other hand, the dressed electron-photon vertex can be 
written in the following from: 



fo(Q)«e 7o + 0(Q), 



(45) 



where e denotes the renormalization of the elementary 
charge due to lattice deformations, such that the effective 
rcnormalized action reads 



<3)*p7o*q- 



(46) 



p Q 



The effect of the lattice defects on the optical con- 
ductivity can be calculated from the dressed polarization 
shown in Fig. [S] The dressing effect of the lattice defects 
is taken into account by replacing bare Greens functions 
G in Fig. S] by the full propagators G defined in Eq. (|4"3")l 
and bare vertices 7 by the dressed ones r o from Eq. 
Algebraically we obtain 

U(K) = J Tr{T G(P)T G(K + P)} 
p 

/z ' PoiPo + k o) - a2 P ■ (P + k ) 



2c: 



[pi + a 2 p 2 ][(p 



A) 2 



a 2 (k + p) 



(47) 



where a = Z± /Zq . The integration can be performed 
with some effort, but we restrict our attempts to a more 
direct task, i.e. we calculate only the modified optical 
conductivity a analogously to the optical conductivity of 
ideal graphene as we did before. We take the derivative 
with respect to k under the integral and employ residue 
theorem in order to integrate out loop frequency p . At 
the end of the calculation we arrive at 



e 2 



Z <7 . 



(48) 



Surprisingly, apart from the vertex renormalization only 
renormalization of the frequency contributes to the mod- 
ificated conductivity. Therefore our task reduces to the 
calculation of the renormalization factors Z and e which 
is performed perturbativelly below. 
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= 1- - — ( E )■*— ' L 




where Tj^ = T 2 = 7 ■ q and 



Q-P E 



FIG. 6: Dressed fermionic Greens function and leading order 
diagram of the fermionic self-energy. Dashed lines denote the 
(AA) -correlator and white circles the ^Avf-vertices. 



b) 



FIG. 7: Second order self-energy corrections at zero external 
momenta. 





F{q) = ¥T^- 



(51) 



denotes the momentum dependent part of the 
(AA)— correlator defined in Eq. ([25)1 . The diagram 
of the next order in g depicted in Fig. [7^, reads 

Z ( i\p ) = ff/ F(q) J F(p)T 1 G (p ,q)T 2 

q p 

xG (p , q + p)T 2 G (p 07 q)r i; 

where T x = 7 ■ q and T 2 = 7 • (2q + p), while the diagram 
shown in Fig. 03 writes 

4 2) (fo) = f|/ *X<?) / F(p)T 1 G (p ,q)T 2 



V. CALCULATION OF RENORMALIZATION 
FACTORS 

In order to set up perturbative calculations we have to 
average over the lattice deformations. There are two pos- 
sible ways to implement such averaging: The replica trick 
and supersymmetry approach. The calculation below is 
based on the replica-trick approach. 

According to the replica trick we introduce N copies 
of fermions \& a , a = 1,2, . . . N with the same action 



S = - 



A P Vf +Q r(P + Q,P,Q)V c c 



Q p 



P+Ql0^Q- 



(49) 



Q p 



Then we will calculate the diagrams describing renormal- 
ization of the fermionic propagator and electron-photon 
vertex function and perform limit N — > at the end 
of the calculation. As the result of this procedure all 
contributions containing factors proportional to any pos- 
itive power of N will vanish. These include for instance 
contributions arising from diagrams containing closed 
fermionic loops. 

The diagrams of the fermionic self-energy to the order 
1 in replica indices and order g, in lattice deformations 
strength are shown in Figs. [6] and [7] Correspondingly, the 
same order diagrams of vertex corrections are depicted in 
Figs. [5] and [SI Retaining only frequency dependence (p = 
0) in the analytical expressions for this contributions we 
obtain for the leading self-energy contribution (Fig. [6]) 



S (2) (Po) 



r 1 G ( Po , 9 )r 2 F( (7 ), 



(50) 



x Go (Po . <? + P) r 3 G (po , p) r 4 , 

with the vertices T x = 7 • q, T 2 — 7 • (2q + p), T 3 = 
7 • (2p + q) and T 4 = 7 • p. Eventually we obtain for 
the contributions to the wave-function renormalization 
factor 



d 
dPo 



S (2) (Po) 



9Po 1 iP0) 

— s (4 V» ) 



I'cr 



P =0 



-9e 2 , 
15 



e 2 . 



where we define 



e=#log^, 



(52) 
(53) 
(54) 

(55) 



with A denoting some upper momentum cutoff. There- 
fore we obtain for the wave-function renormalization to 
the second order in g% 



Z Q - 1 a 1 



fe 2 + C(e 3 ). 



(56) 



Now we look at the renormalization of the electron- 
photon vertex function. Diagrammatically the leading 
order correction is given by the second term on the right 
hand side of the diagram shown in Fig. [8] According to 
Eq. (|48p the main corrections to the conductivity arise 
from the momentum independent part of the vertex func- 
tion. We obtain for the vertices T l = T 2 = 7 • q. The 
algebraic expression for the diagram depicted in Fig. [8] is 
given by 



f«(0) 



9l 



r 1 G (q)r G (q)r 2 F( 9 ), (57) 
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FIG. 8: Leading order perturbative contribution to the ampu- 
tated electron-phonon coupling vertex. Dashed lines denote 
the (AA)-correlator and white circles the ^A^-vertices. 




FIG. 9: Second order vertex correction diagrams at zero ex- 
ternal momenta. The diagrams in the second row should be 
counted twice due to the mirror symmetry. 

which yields after the evaluation 

fW(0) = e 7o , (58) 

where e from Eq. (1551) is introduced and acquires the 
meaning of the leading order elementary charge renor- 
malization. 

Second order vertex corrections can be calculated from 
the diagrams shown in Fig. [§] Due to the mirror sym- 
metry diagrams c) and d) depicted in the second row 
in Fig. [5] should be counted twice. Let us first consider 
the contribution that arises from the diagram depicted in 
Fig. [9^ with two parallel ladder rungs. At zero external 
momenta we obtain for the vertices T 1 = T 4 = 7 • q and 
r 2 = T 3 = 7 • (2q + p) . We obtain for the correction 

ff(0) = F(q) J F(p)r l G (q)T 2 G Q (p + q ) 

q p 

xT G {p + 0)^0(0)^, (59) 
which, after performing integrations, yields 

ff>(0) = 2e 2 7o . (60) 

For diagrams depicted in Figs. [9p, c and d we proceed 
similarly. In the case of diagram b) we have the following 
expressions for the vertices: T 1 — 7 • q, T 2 — 7 ■ (2q + p), 
T 3 = 7 - (2p + q) and T 4 = 7 -p. Therefore the expression 



for this correction reads 
ff^O) = F(q) J F(p)T 1 G () (q)r 2 G (p + q) 

q p 

xT G (p + q)T 3 G (p)T 4 , (61) 
with the result 

f^ 2) (0) = ^ 2 7o . (62) 

For diagram c) we have the following vertices: T 1 =7-9, 
T 2 = 7- (2q+p), T 3 = 7 - (2p + q) and T 4 = 7 -p, whereas 
the expression for the correction reads 

if } (0) - f|y F(q) J F(p) r 1 G (g)r G (g)r 2 

q p 

xG Q (p + q)T 3 G (p)r 4 , (63) 
which yields the result 

ff (0) = ^e 2 7o . (64) 

Finally, diagram d) from Fig. [5] can be written alge- 
braically as follows: 

rf(0) = f|y F(q) J F(p) r 1 G (g)r G (g)r 2 

q p 

xG (p + q)r 3 G (q)r 4 , (65) 

where the vertices T^.^ are given by T 1 — T 4 = 7 ■ q and 
r 2 = T 3 = 7 • (2q + p). Evaluation of Eq. (JB3J) yields 

fi 2) (0) = ^e 2 7o . (66) 

Hence, the second order contribution to the vertex 
function becomes 

f (2) (0) = if ] (0) + f 2 2) (0) + 2ff ) (0) + 2lf ] (0) , (67) 

such that the dressed vertex function can be written as 
a series in e: 

f o = e 7o w e (l + e + ^e 2 + 0(e 3 )^j 7o . (68) 

This expression reproduces exactly the result which we 
have obtained for the dressed vertex function in Eq. (|56l) . 
Therefore we obtain from Eqs. (|48|) . (|68|) and (|56| for the 
modified conductivity 

er« {1 + 0(e 3 )) a 0l (69) 

i.e. the leading correction is of the order <7q. However, 
we can show to every order in perturbative expansion 
that corrections arising from the propagator renormaliza- 
tion are exactly canceled by their counterparts departing 
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from the electron-photon vertex renormalization. Con- 
sider the definition of the quasi-particle weight Z^ 1 given 
in Eq. (J44j . Since the lattice deformations are static, all 
propagators inside the diagram depend only on the ex- 
ternal Matsubara-frequency which thus becomes an in- 
dependent parameter. Hence the derivative with respect 
to the Matsubara frequency should be applied to every 
propagator. Taking such a derivative of an average free 
propagator 

(^o^%=f 9 °y ( ?! + ^ , (70) 



gg + (q + kf 

k 

at zero external momentum and frequency we obtain 
d 



0q { 



■(G {q ,q + k)) k 



-<G (0,fc)7cA(0,fc)) fc .(71) 



Q=0 



Equation (I7T1) suggests that the expressions under the 
integrals must be equal up to an irrelevant constant. 
Therefore each derivative of the free propagator with re- 
spect to the external frequency generates upon sending 
external momenta and frequency to zero a bare electron- 
photon vertex. An irreducible n— th order diagram of the 
electronic self-energy contains 2n — 1 electronic propaga- 
tors. Therefore by applying a derivative with respect to 
the frequency q to such a diagram, 2n — 1 irreducible 
corrections to the electron-photon interaction vertex are 
generated. This mimics term by term a perturbative se- 
ries for the dressed vertex. This can be seen very clearly 
if we look at self-energy diagrams depicted in Fig. [7] 
Replacing electronic propagators successively by a bare 
electron-phonon interaction vertex and putting external 
momenta to zero we reproduce exactly vertex-correction 
diagrams shown in Fig. |H1 We therefore can link each n- 
th (n 1) term in the perturbative series of self-energy 
to the corresponding vertex function correction: 



d 



(Q) 



-f{, ,l) (0,0,0). 



(72) 



2=o 



Summing over all n and subtracting 7 on both sides we 
then can assemble all contributions arriving at 



= f (0,0,0). 



(73) 



On the right-hand side we have the charge renormaliza- 
tion e7 defined in Eq. ([45]) , while the left-hand side rep- 
resents the wave- function renormalization factor Zj~ 1 7 
due to Eq. (|43")) . Therefore Eq. (1731 postulates the equal- 
ity 



e, 



which leads to the exact result for the modified conduc- 
tivity: 



Importantly Eq. (|73[) is obtained without special empha- 
size on a disorder type and is not restricted to the con- 
sidered type. The only requirement we need is that the 
corresponding term should not violate the chiral sym- 
metry of the pure graphene Hamiltonian and must be a 
quenched disorder. 



VI. CONCLUSIONS 

In the present paper we address the question of the ef- 
fect which random deformations may have on the trans- 
port in graphene. The common believe is that surface 
corrugations in graphene influence its electronic trans- 
port properties, mainly the optical conductivity. It is 
possible to describe deformations in graphene by a gauge 
field that couples to the fermions living on the two di- 
mensional sheet. We have performed perturbative calcu- 
lations of the corrections due to lattice deformations to 
the optical conductivity. Our results contrast the sugges- 
tions made in Refs. dUll where a substantional effect 
of the defects on the conductivity is proposed. We have 
found that the minimal conductivity is robust with re- 
spect to the surface corrugations. 
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Appendix 

Below be evaluate the irreducible polarization of ideal 
graphene starting with Eq. (|40l) . Upon performing the 
trace over the pseudo-spin space we arrive at 



n(K) 



2 Po (p + k ) - P ■ (P + K) 



where P ■ K = p k 
parametrization 

1 

AB = 



P 2 (P + K) 



p ■ k. Employing the Feynman 



dx 



[xA + (1 -x)B] 2 ' 



and shifting P — > P—xK we symmetrize the denominator 
with respect to P. Therefore odd powers of P appearing 
in the numerator may be dropped. We arrive at 



(74) 



IL(K) 



dx 



2pl-P 2 + x{\-x)[K 2 -2k 2 } 
[P 2 + x(l - x)K 2 } 2 



9 



Exploiting the rotational invariance we replace p 2 , = which yields after integrating out x the result of Eq. (|41l) 
P 2 /3 and use formulas of dimensional regularization: 

d d k 1 1 r(n-f) 



/ 



(27r) d (k 2 + A)" (4tt)<V 2 r(n)A»- d / 2 ' 
d d fc k 2 1 d T(n - f - 1) 



(27r) d (fc 2 + A)™ ' ' (47r) d / 2 2 T(n)A n ~ d / 2 - 
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